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Abstract 

Observations of diffuse Galactic 7-ray spectrum by the EGRET instrument reveal an 
excess above ~ 1 GeV over the expected 7-ray spectrum calculated under the assump- 
tion that the locally observed cosmic -ray (CR) spectra represent the galactic CR spectra. 
Assuming that Galactic CRs of energy below ~ 100 TeV are accelerated by supernova 
remnant (SNR) shock waves and that the shock compression ratio is SNR age dependent, 
the average source injection spectra from an ensemble of SNRs is calculated both in the 
inner (330° < I < 30°) and outer (30° < I < 330°) regions of the galaxy. The calculation 
considers the SNR age distribution in the galaxy. Injecting these spectra in the galaxy and 
using a 3-D convection-diffusion equation, the CR electrons and protons spectra in the two 
galactic regions are obtained and their spectra in the galactic disk are found to be flatter than 
the observed CR spectra. The diffuse gamma-ray spectrum produced by the interaction of 
these galactic CRs with the ISM and ISRFs is compared with the experimental data in both 
the galactic regions. Furthermore, the steepening of the observed local CR spectra from the 
galactic disk CR spectra are discussed by propagating local CRs having a source spectrum 
derived using local SNR age distribution (SNRs located within l.bkpc from the Sun), for a 
diffusion coefficient Dq ~ 0.3 x 10 2 cm 2 s~ 1 in the local region which is much less than 
the typical value in the galaxy Do ~ (1 — 10) x 10 28 cm?s~ l . The results obtained in this 
paper support the SNR origin of galactic CRs. 
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1 Introduction 



The spectrum of radiation produced by the interactions of galactic CRs with the 
interstellar medium (ISM) and the interstellar radiation field (ISRF) gives informa- 
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tion about their propagation in the galaxy, their sources distribution and their ini- 
tial source spectrum. Results from observations made with the Energetic Gamma 
Ray Experiment Telescope (EGRET) on the Compton Gamma-Ray Observatory 
have shown that above ~ 1 GeV, the diffuse Galactic 7-ray spectrum exceeds sig- 
nificantly the spectrum calculated using the locally measured proton and electron 
spectra [1]. Explaining the excess with the 7r°-decay gamma-ray alone requires a 
proton spectrum with index ~ 2.4—2.5 [2] which is much flatter than ~ 2.7, the 
locally measured value [3]. Many authors have given different views to explain the 
excess (see e.g. [4-9]), but, still the problem is not fully understood. 

Although supernova explosions are considered to be the main sources of CRs with 
energies up to ~ 10 14 eV in the galaxy (see e.g. [10]), most of the SNR statistics 
(except the radio spectral index distribution and the source density radial distribu- 
tion) have not been used extensively in studying this puzzle of "GeV excess" and 
the observed cosmic-ray spectra. In most of the earlier studies (e.g. [4,5]), the CR 
source injection spectrum was taken either from the radio observations (taken at 
a few GH Z which corresponds to an electron energy ~ 5GeV) or as an input pa- 
rameter which was then optimized to the observed 7-ray data. In this paper, the 
injection spectrum is first calculated independent of the 7-ray data by considering 
shock acceleration mechanisms and then injected into the galaxy and see whether 
the excess can be explained. This is more realistic compared to the earlier studies. 

In spite of the fact that ~ 77% of the total 231 galactic SNRs detected in radio 
are shell(S)-type, statistical studies have found that there is no firm relationship 
between the radio spectral index and the age of the remnant [11,12] while ideally 
shock acceleration theory predicts a correlation between the particle spectral index 
q(t) (which relates with the radio index) and the shock velocity (which in turn 
relates with the SNR age) as [13] 

aft) = = mt)2 (1) 
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where s(t) = , is the shock compression ratio, M(t) = u(t)/c s is the Mach 

number of the shock defined as the ratio of the shock velocity u(t) at an age t to 
the sound speed of the unshocked intersteller medium c s and 

u(t) = m , t < tsed, (la) 
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t S ed = (3M ej / AirfXTriHriHul) 1 ^ 3 denotes the time which marks the start of the Sedov 
phase where M ej is the ejecta mass, jj, = 1.4 is the mean atomic weight of the 
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ISM, m H is the hydrogen mass, n H is the mean ISM hydrogen density and u = 
^J2Esn I 'M e j is the initial shock velocity with E$n repesenting the energy released 
in SN explosion. t ra d = 2.7 x lO^E^nJj ' 52 yr is the start of the radiative phase 
which also represents the end of the Sedov phase where E 51 is the SN explosion 
energy in units of 10 51 ergs [14]. The non-correlation between the radio index and 
the SNR age requires a very detailed study and will not be discussed here, but one 
possibility may be due to the presence of more than one electron population, as is 
well known in the case of Crab nebula. 

In section 2 an estimate of the SNR age distribution in the galaxy is given. This 
distribution is then used in section 3 for calculating the particle source spectra in 
the inner (330° < I < 30°) and outer (30° < I < 330°) regions of the galaxy. 
The resulting source spectra are further used to obtain the respective galactic CR 
spectra. In section 4, using the galactic CR spectra the diffuse gamma-ray flux is 
calculated. In section 5, the calculation of the local CR spectra is given, and in 
section 6, the over-all results and the limitations of the model are discussed briefly. 
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Fig. 1. Normalised distribution of 77 S-type SNRs with known ages in the galaxy. The 
thick curve is the fitted gaussian (equation 2) with a = 0.266 ±0.045, ai = 4. 147 ±0.104 
and a 2 = 0.831 ±0.113 



2 SNR age distribution 



Table 1 in Jian et al. 2005 [12] had given a list of several parameters of 23 1 all-type 
SNRs, out of which only 77 S-type SNRs are with known ages which have been 
estimated by different methods as indicated in the list. The normalized distribution 
of these 77 S-type SNRs (known ages) is shown in Fig.l. The distribution can be 
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fitted with the gaussian 




x — d\ 
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n(x) = a exp —0.5 



(2) 



where x = log w t and t is the SNR age in years. The fitted parameters together 
with their statistical errors, reduced chisquare (x 2 /no. of degrees of freedom) and 
fit probability (Prob) are given in Table 1. In the present study, it will be assumed 
that the age disrtribution in Fig.l represents a galactic averaged S-type SNR age 
distribution in the galaxy. 

It can be argued that the age distribution shown in Fig.l may not represent the 
true age distribution in the galaxy because of selection biases in radio observations 
which arises mainly due to the limited sensitivity, field of view and angular res- 
olution of the instruments. Among these, only the limitation due to instrumental 
sensitivity is likely to affect the SNR age distribution strongly since it is directly 
related to the brightness of the source. Other effects are only likely to reduce the 
sample size as a whole irrespective of the SNR age. But, since the surface bright- 
ness decreases with SNR age, the sensitivity biasing will have their effects mostly 
on old SNRs. In fact, assuming adiabatic phase the surface brightness estimated 
from an SNR of age ~ 10 5 yrs is fainter than the detection limit in studies of 
Galactic SNRs distribution [15,16]. For the present study, we are mainly interested 
in SNRs with ages between approximately 10 3 yrs and 3 x 10 4 yrs (the start and 
end of the Sedov phase respectively, see next section) since particle acceleration 
is assumed to occur mainly in the Sedov phase of the SNR evolution. So, it will 
be assumed that the effects of detection biases due to the finite sensitivity of radio 
telescopes are quiet small in the present study. One more approximation concerns 
the confinement time of CRs in SNRs. Detailed treatment of particle acceleration 
in supernova shock waves involves the confinement of accelerated particles in the 
shock region until the shock becomes weak. Although the most energetic particles 
start escaping from the acceleration site already at the start of the Sedov phase, the 
main fraction of particles remain confined in SNRs for approximately (10 4 — 10 5 ) 
yrs [17]. But, for simplicity, it is assumed here that after acceleration (which starts 
at the onset of the Sedov phase) the particles start escaping rapidly from the shock 
region into the ISM. Therefore, from now on it will be assumed that the age distri- 
bution shown in Fig.l also represents the distribution of CR acceleration times in 
the galaxy. 

One interesting point that can also be noted is that the total spectrum from an en- 
semble of SNRs is dominated mainly by particles liberated from SNRs which are 
near the Sedov phase, except at the very low energy end of the spectrum (near parti- 
cle injection energy ~ 1 KeV). It is because of the high acceleration efficiencies of 
young SNRs which can give rise to flatter particle spectra in contrast to low accel- 
eration efficiencies of old SNRs which give them a comparatively steeper source 
spectra. So, it can be concluded that almost the whole range of the total particle 
spectrum depends strongly on t sed and very weakly on t rad . A typical hydrogen 
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Table 1 

Fit parameters for the normalized age distribution of 77 S-type SNRs with known ages 
shown in Fig.l : 



Entries 


11 


X 2 /ndf 


3.258/4 


Prob 


0.516 


a 


0.266±0.045 


ai 


4.147±0.104 


a 2 


0.831±0.113 



density of tih ~ 0.5cm" 3 gives t sed ~ 1.3 x 10 3 yrs and t rad ~ 3.8 x 10 4 yrs, both 
of which lie well within the detection limit of radio telescopes. Therefore, the ef- 
fects of detection biases on the final source spectrum are definitely small. However, 
for a hydrogen density much lower than the typical galactic averaged value, say 
n H ~ 0.01cm -3 , t sed ~ 4.7 x 10 3 and t rad ~3x 10 5 . Here also, although t rad lies 
below the detection limit, the effect of the SNR selection biases on the total source 
spectrum will be considerably small as long as t sed lies well within the sensitivity 
of radio telescopes since young SNRs whose ages are around t sed dominate the fi- 
nal spectrum. But, such an extreme condition is far away from the typical galactic 
conditions (n H ~ 0.5cm -3 ) which is considered in this work. 



3 Galactic cosmic-ray spectrum 

3.1 Source CR spectrum 

Considering the similarity of the power supplied by a SNe (10 51 erg/30yrs ~ 
10 42 erg/s) to the power required to maintain the CR energy density in the galaxy 
(~ 10 41 erg/s), it is widely beleived that the majority of CRs upto ~ lOOTeV are 
accelerated in SNRs. There is a wealth of literatures available which discusses the 
theory of diffusive shock acceleration both in a simple planar shocks model (see 
e.g. [18-20]) and in more detailed models that consider the shock goemetry as well 
as the non-linear CR backreaction effects (see e.g. [21-23]). In the present study, 
the non-linear effects are neglected and the much simpler planar shocks model will 
be adopted. 

Acceleration of cosmic rays in plane, steady shocks in which the cosmic rays do 
not influence the shock structure produces a power law spectrum of the form 

f(p,t)=q(t)p-^ !* dp gtp ^p™- 1 (3) 
Jo 
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where f(p,t) is the isotropic accelerated particle distribution function with mo- 
mentum p arising from an injected particle distribution g(p ,t) with momentum p 
at time t. For constant injection of monoenergetic particles of momentum p with 
density n , 

g(p',t) = -^S(p-p ) (4) 
which gives the accelerated particle spectrum at time t as 

fM = £^q(t)(p/po)- q V (5) 

Assuming that the SNR age distribution shown in Fig.l exists in the galaxy at 
any point of time and injects CRs continuously in the galaxy, the particle source 
spectrum from an ensemble of SNRs is calculated by multiplying equation (5) by 
the SNR age distribution function (equation 2) and integrating from t = t sed (the 
start of the Sedov phase) to tf = t rad (the end of the Sedov phase) as 

f(p) = ^3 f Q(t)(p/Por q{t) Hx)/t]dt (6) 
Anpl Jto 

For calculating the actual particle source spectrum in the galaxy, the averaged ISM 
hydrogen atom densities in the inner and outer regions of the galaxy are calcu- 
lated as n H = n H (R)dR/ (R 2 -Ri) where R is the galactic radius, n H (R) = 
fa n H (R, z)dz/h and n H (R,z) = n H i(R,z) + 2n H2 (R,z). Here, (Ri,R 2 ) are 
taken as (2.25, 8.5)kpc and (5, 16)kpc for calculations in the inner and outer galaxy 
respectively and z represents the galactic height. The atomic hydrogen density dis- 
tribution tihi(R, z) is taken from [24] and the molecular hydrogen density n H , 2 (R, z) 
is from [25]. Different authors had given slightly different values for the CO-to-H 2 
conversion factor X = N(H 2 )/W(CO) where N(H 2 ) is the molecular hydrogen 
column density and W(CO) is the velocity integrated CO intensity. X = 3.6 x 10 20 
cm' 2 (K km s' 1 )' 1 by Sanders et al.1984 [26], X = (2.8 ± 0.7) x 10 20 cm' 2 (K 
km s' 1 )' 1 by Bloemen et al. 1986 [27] and X = 1.9 x 10 20 cm' 2 (K km s' 1 )' 1 
by Strong & Mattox 1996 [28]. Here, the value given by Strong & Mattox(1996) 
is adopted. Such a distribution of n HI (R, z) and n H , 2 (R, z) gives a mean hydrogen 
atom density of n H ~ 0.5(0.3) cm' 3 in the inner(outer) galaxy for z < 700pc. 
Correspondingly, this gives a value of t sed = 1.3 x 10 3 yr & t rad = 3.8 x 10 4 yr in 
the inner region and t sed = 1.5 x 10 3 yr & t rad = 5 x 10 4 yr in the outer galaxy. Us- 
ing these values and considering a total energy released of E S n = 10 51 ergs in SN 
explosion together with a total ejecta mass of M ej = 8M© where M@ is the solar 
mass, the CR proton and electron source spectra can be calculated from equation 
(6) for the normalized age distribution shown in Fig. 1 . The resulting energy spectra 
are plotted in Figs. 2&3 for the inner (330° < / < 30°) and outer (30° < I < 330°) 
regions of the galaxy respectively for the case when particles are injected at energy 
1 KeV with number density no = 10~ 3 cm~ 3 . Detailed inspection of Figs. 2&3 
reveals the following: 

(i) For protons, the energy spectrum follow a broken power-law 1(E) oc E' Fp with 
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Fig. 2. CR source spectra in the inner region (330° < I < 30°) of the galaxy calculated 
assuming E$n = 10 51 ergs, M e j= 8M(g, Tin = 0.5cm~ 3 , uq = 1CT 3 cm~ 3 , injected 
particle energy of 1 KeV and an upper energy cut-off of 100 TeV. Proton source spectrum 
(solid line): spectral index T p ~ 2. 12 ±0.01 for E » M p and 1.66 ±0.01 for E « M p . 
Electron source spectrum (dashed line): T e ~ 2.126 ± 0.004 
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Fig. 3. CR source spectra in the outer region (30° < / < 330°) of the galaxy. All the 
parameters are same as in Fig. 2 except fiu = 0.3cm -3 . Proton source spectrum (solid 
line): T p ~ 2.11 ± 0.01 for E » M p and 1.66 ± 0.01 for E « M p . Electron source 
spectrum (dotted line): T e ~ 2.124 ± 0.004 

a break at ~ M p , where E is the kinetic energy and M p is the proton rest mass 
energy. The spectral index T p ~ 2.12 ± 0.01 for E » M p and 1.66 ± 0.01 for 
E « M p for the inner galaxy and 2.11 ± 0.01 for E » M p and 1.66 ± 0.01 for 
E « M p for the outer galaxy. 

(ii) For electrons with E » M e where M e is the electron rest mass energy, the 
spectrum has a power law index Y e ~ 2.126 ± 0.004 for the inner region and 
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2.124 ± 0.004 for the outer region. 
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Fig. 4. Plot of q(t) vs logwt for the inner galaxy, where t is the SNR age in years. The solid 
dots (•) marked as 2,5,9&12 are the values of q(t sed ) = (4.011,4.028,4.051,4.068) for 
M e j = (2, 5, 9, 12)M(g respectively. The sedov phase ends at t rad = 3.8 x 10 4 yr which is 
not shown in the plot. The dashed curve is the particle momentum index q(t) obtained from 
equations (1) &(lb) for SN explosion energy E$n = W 51 ergs. 



These spectra are almost consistent for a wide range of ejecta masses. For exam- 
ple: In the inner region, the spectral index changes from (2.091—2.136) for protons 
with E » M p and from (2.102-2.146) for electrons, i.e. by a value of ~ 0.044 
for ejecta masses (2 — 12) M® for the same explosion energy of 10 51 ergs. This 
mass range covers a large number of both the Type I and Type II supernovae. This 
consistency can be understand as follows: When the total spectrum from an en- 
semble of SNRs is calculated from equation (6), because of their high accelera- 
tion efficiencies which gives rise to a comparatively flatter spectra, those which are 
near the Sedov phase (t sed ) dominate the final spectrum except at the very low en- 
ergy end of the spectrum (near particle injection energy ~ 1 KeV) as mentioned in 
the last section. Fig. 4 shows the q(t) vs t plot for the sedov-phase obtained from 
equations (1) &(lb). Since u oc M~^ 2 and t sed oc (Mej/wq) 1 / 3 , it is found that 

t sed increases with M ej as t sed oc M^ 6 . Also shown in the figure are the values 
of q{t sed ) = (4.011,4.028,4.051,4.068) for M ej = (2, 5, 9, 12)M @ respectively 
which changes by a small value of ~ 0.057 for the specified mass range. This weak 
dependence of q(t sed ) on M e j is expected since q(t sed ) = 8E$n /(2Esn — c sMy ) 
from equations (1) &(lb). This slow variation of q(t sed ) is responsible for the small 
changes ~ 0.044 in the spectral index of the final particle spectrum obtained from 
equation (6) which is the superposition of individual spectra from various SNRs 
with ages (t sed < t < t rad ). 

It should be mentioned at this point of discussion that the mean hydrogen atom 
density obtained here may be less by a factor of upto ~ 2 than that used for cal- 
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culating the diffuse gamma-rays from a particular direction in the galaxy because 
of the dominant emission of the radiation fluxes from a region very close to the 
galactic plane, say within a region of z < lOOpc [6]. Thus, for calculating the dif- 
fuse gamma-ray fluxes later in section 4, an averaged atomic hydrogen density of 
~ l(0.6)cm~ 3 will be adopted in the inner(outer) galaxy. 



3.2 Propagation ofCR protons in the galaxy 

In the framework of a stationary 3D convection-diffusion model, the propagation 
of cosmic ray protons in the galaxy can be represented by the equation 

V ■ {B(r, E)VN p (r, E) - V(r)JV p (r, E)} 

+ we { (™ B - ^ E) ) N ^ E) } - 3£f + Q ' [T ' E) = (7) 

where E is the proton kinetic energy in GeV, N p (r, E) is the differential number 
density , D(r,E) = D(E) = D (l + E/E ) a [29] is the diffusion coefficient 
which is assumed to be independent of r with a = 0.6, E = 3 GeV and D = 
(1 - 10) x 10 28 cm 2 s~ 1 , Q p (r, E) is the proton production rate, i.e. Q p (r, E)d 3 rdE 
is the number of protons produced by the source in a volume element d 3 r in the 
energy range (E, E + dE) per unit time, V(r) = 3V zz is the convection velocity 
in the direction perpendicular to the galactic plane with V = ISkms^kpc^ 1 and 
P p (r, E) = —Bp(E) is the proton energy loss rate due to ionization and Coulomb 
interactions. Following Mannheim and Schlickeiser 1994 [30], if E is in GeV and 
proton mass energy, m p = 0.938 GeV, then 

BJE) w 1.82 x 10" lf W (—] (l + 0.85^-) GeV s' 1 (8) 

\m p J \ n H J 

where n HII is the mean density of ionised gas in cm~ 3 , [—N p (E : r)/r pp (E,r)], 
t pp = E/(dE/dt) v and (dE/dt) n = 8 x 10 _16 n H E in GeV s" 1 are the catas- 
trophic loss term, energy loss time scale and energy loss rate respectively due to 
pion production processes. With these substitutions, equation (7) can be simplified 
in rectangular coordinates (x,y,z) as 

( d 2 N d 2 N d 2 N 1 d d 

D ^ {k + W + St) - m ^ zN ' } + ^ {{V ° E - ME)) N ' ] 

N 

-^ = -Q P (r,E) (9) 
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The Green function of equation (9), i.e. the function G(x, x , y, y\ z, z , E, E') sat- 
isfying the equation 



D(E) 



d 2 G d 2 G d 2 G 



| - SV ^(zG) + A { {y Q E - (3 P {E)) G} 



d 2 x d 2 y d 2 z 

-— = -5(x - x')S(y - y')S(z - z')5(E - E') (10) 



'VP 



is found to be 



G(x,x,y,y,z,z,E,E) = 87r3/2pjE) ^ 
rE' D{u)A 2 {u)du 



E ' D(u)du 



f 
Je 



P p (u) 



P (E) [Je P p {u) 



rE du 



xexp 



- -1/2 


f 1 


exp 


A pp J e 






-((x -x 


) 2 + (y 



xexp 



a rE 1 D(u)du 
4 JE P v ( u ) 



. rE' D(u)Al(u)du 
^JE Al(E)P p (u) 



where P P (E) = V E + B P {E) and 



A P (E) = exp 



3Vb 



du 



V u + B p (u)_ 



(ID 



x 10~ lf W 



Then, the general solution of equation (9) is obtained as 



/■oo , roc , roo , roo , , , , , 

N p (x,y,z,E)— / dx I dy dz dE Q p (x ,y , z , E ) 

J — oc J— oo J— oo J E 

xG(x,x',y,y',z,z,E,Ef) (12) 

Now, considering a source of the form Q(x ,y , z , E') = q(E')S(x — x )8(y' — 
yo)5(z — z ), the differential proton number density at (x,y,z) due to a single 
source located at (x , y , z ) is obtained as 
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- E ' D{u)du 



N p (x, x , y, y , z, z , E) = Je [Je P p (u) 

r& D(u)A 2 (u)du 



Je 



P p (u) 
xexp 



1-1/2 




exp 


A pp J e 







-((x -x) 2 + (y -y) 2 ) 



xexp 



a rE' D(u)du 
4 j. 

-(3V^ jf- Ap{u)d " 



IE 1 A p (E)P p (u) 



- Z + Z 



. r E' D(u)Al(u)du 
^JE A$(E)P p (u) 



(13) 



This gives the number density at (x — 0, y — 0, z) as 



f 
Je 



exp 



r E ' D(u)du 
P P {u) . 
E du 



-i 



PP 'e' P p (u) 



exp 



-(4 + y 2 o) 




exp 


, rE' D(u)du 
l^JE Pp (u) _ 



P p (u) 

(ZVqZq J§ Ap( E)p d p U (u) - ZQ + z) 



21 



. r-E' D{u)Al(u)du 
4JE A*(E)P p (u) 



(14) 



which in cylindrical co-ordinates (r , 9 , z ) becomes 



X 



-1/2 



exp 



f E du 
'e' P p {u) 



-i 



A 



exp 



—r'r 



a rE' D{u)du 
^JE Pp (u) 



exp 



- {ZVqZq J e ' a p {e)p p {u) ~ z o + z 

a rE' D{u)Al{u)du 
^JE Al(E)P p (u) 



21 



(15) 



For a uniform source distribution in the galactic disk with density F(r , 9 , z ) = 
kS(zo), the number density N(E) at (0, 0, z) due to all the sources located in the 
disk region with radial distances between r 1 and r 2 from (0, 0, 0) can be calculated 
by integrating equation (15) over the galactic region bounded by r\ and r 2 as 



11 



N p (z,E) 



kA p {E) 



2^P P (E) Je 

rE du 



dE'q(E 



xexp 



A 



pp 



f 

Je' 



P p (u)\ 




exp 



xexp 



. rE 1 D(u)Al(u)du 
^JE A^E)P p (u) 



(16) 



where 



rE 

Je 



E D(u)du 



This equation represents the variation of the proton density with the galactic height. 
For a CR source distribution in infinite space, by setting z = equation (16) rep- 
resents the CR proton density at any point in the galactic disk due to sources lo- 
cated within the galactic disk region bounded by r x and r 2 from the point. Taking 
r 2 — > oo, equation (16) gives the proton density at (0, 0, z) due to all the sources 
located in the disk beyond a distance r x from (0, 0, 0) as 



N p (z,E) 



xexp 



A 



kA p (E) 
2^P p (E) Je 

E du 



pp 



f 

Je' 



dE'q(E') 



exp 



E ' D(u)A 2 Ju)du 



E 



exp 



P p (u) 



-1/2 



. r E' D{u)Al(u)du 
4 Je AUE)P p {u) 



(17) 



In the case when r\ — > 0, equation (17) gives the number density at any galactic 
height z due to all the sources present in the galactic disk as 



N p (z,E) = 



kA p {E) 
2^P P {E) Je 



xexp 



A 



dE'q(E' 



du 



e' D(u)Al(u)du 



2 ( 
P p (u) 



-1/2 



PP V P P (u) 



exp 



. rE' D{u)Al(u)du 
^JE Al{E)P v {u) . 



(18) 



3.3 Propagation ofCR electrons in the galaxy 



The form of the transport equation for CR electrons in the galaxy differs from that 
of the protons only in the absence of the catastrophic loss term as 
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V • {D(r, E)VN e (r, E) - V(r)iV e (r, E)} 

+ ^e { {^T~ E ~ /3e(r ' E) ) Ne{r > E) ] + ge(r ' E) = (19) 

Here, /3 e (r, E) = —B e (E) represents the total electron energy loss rate due to ion- 
ization, bremsstrahlung, inverse compton and synchrotron processes. From Atoyan 
etal. 1995 [29], 

B e (E) = 3.07 x l(T lf W + 1.0 x lO~ 15 n H E + 

1.01 x 10- 16 (w ph + w B )E 2 (20) 

in GeV s" 1 . Here, n H is the mean hydrogen atom density in atoms cm" 3 , w ph = 
wmbr + w NIR + w FIR and w M br, w nir , w F ir, w b are the energy densities of the 
microwave background radiation, difuse NIR/optical radiation, diffuse FIR radia- 
tion and magenetic field respectively in eV cm" 3 . Following the same mathematical 
calculations given in section 3.2 for protons on the assumption of sources uniformly 
distributed in the galactic disk, the electron density N e (z, E) at (0, 0, z) due to all 
the sources located in the galactic region bounded by the radial distances r 1 and r 2 
from (0, 0, 0) is obtained as 



r& D(u)A 2 e (u)du 



PAu) 



-1/2 



exp 



— exp 



exp 



a rE' D(u)Al(u)du 
4 Je Al(E)P e {u) 



where P e (E) = V E + B e (E) and 



(21) 



A e (E) = exp 



SV 



du 



VqU + B e (u) 



E D(u)du 



"' " ' h: PA a ) 



The electron density at (0, 0, z) due to all the sources located beyond a distance r\ 
from (0, 0, 0) in the galactic disk is given by the equation 



N e (z,E) 



kA e {E) 
2^P e (E) Je 



dE'q(E' 



7 ' D{u)A 2 e (u )du 
Pp.(u) 



1-1/2 



xexp 



exp 



a rE' D(u)Al{u)du 
4 JE 



(22) 



Al{E)P e {u) 
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and the number density at any galactic height z due to all the sources present in the 
galactic disk is given by 



N e (z,E) 



kA e {E) 
2^P e (E) Je 



dE'q(E' 



e' D(u)A 2 e (u)du 

- Pe(u) 



1-1/2 



xexp 



a rE' D{u)Al(u)du 
^JE Al(E)P e (u) 



(23) 




Fig. 5. Proton spectra in the inner (solid line) and outer (dashed line) galactic disks (z = 0) 
calculated in the framework of convection-diffusion model using equation (18) for an in- 
jection spectrum given by equation (6). The spectra are normalized to the observations at 
WGeV. The calculation assumes nn = 0.5(0.3)cm~ 3 and nun = 0.05(0.03)cm -3 for 
inner(outer) galaxy respectively. It can be seen that the two spectra almost overlap each 
other when normalised to observations at WGeV. Data points : ■ Ryan (1972) [31], A Al- 
caraz 2000 [32], x Belloti 1999 [33], ▲ Sanuki 2000 [34], V Boezio 1999 [35], o SOKOL 
[36], T JACEE [37] and • MUBEE [38]. 



3.4 Calculated particle spectrum compared with observed spectrum 



In the convection-diffusion model, for CR protons and low energy electrons which 
do not suffer significant radiative losses we can define a characteristic scale length 
z c which is basically the height above the galactic plane at which the charac- 
teristic diffusive time scale ~ z 2 /D(E) equals the convective timescale V G l as 
z c = [2D(E) /(3 + a/2) Vo] [46,47]. The effect of convection becomes visible only 
when z c « z h , the halo boundary which is taken as infinity in our model. z c is 
essentially the boundary between a diffusion dominated region (z < z c ) and a con- 
vection dominated region (z > z c ) and in fact the boundary of an effective diffusion 
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Fig. 6. Electron spectra in the inner (solid line) and outer (dashed line) galactic disks 
(z = 0) in the framework of convection-diffusion model calculated using equation 
(23) normalized to observations at WGeV for an injection spectrum given by equa- 
tion (6). The calculation assumes uh = 0.5(0.3)cm~ 3 , uhii = 0.05(0.03)cto~ 3 , 
B = 10(6)/xG and wnir = 2. 0(0. 5)eV/ cm 3 for the inner(outer) galaxy. For both the 
regions wmbr = 0.25eV/ cm 3 and wfir = 0.2eV/cm s . Data points : + Alcaraz 2000 
[39], o Boezio 2000 [40], x Du Vernois 2001 [41], A Grimani 2002 [42], ■ Nishimura 
1980 [43], ▲ Nishimura 2001 [44], □ Barwick 1998 [45]. 

halo beyond which the CR density decreases strongly. It plays the same role as the 
diffusion halo in the pure diffusion models except that it is energy dependent. For 
the diffusion dominated regime i.e. z « z c , say at z = 0, N(E) oc E~( v+a l 2 \ 
This spectrum is flatter than the one expected from pure diffusion models which 
can be obtained here by setting Vq — > in our model. For this, z c — > oo from the 
above relation which implies z c ~ z h and the effect of convection becomes negli- 
gible. Then, our model approaches that of the pure diffusion models which gives 
N(E) oc E~( v+a \ Similar characteristic scale length can also be defined for high 
energy electrons which suffer significant radiative losses (with rates given by — bE 2 
where b is a constant) as za = [AD(E)/((\ — a)bE)\ [46]. But here, the compet- 
itive timescales which define the scale length are the diffusive timescale and the 
radiative losses timescale. For such high energy electrons, the electron distribution 
approaches that of the pure diffusion model with radiative losses and is given by 
N{E) oc £-(rW2+i/2) for z <K Zd _ 

The spectrum of protons and electrons calculated using equations (18)&(23) re- 
spectively for the inner(solid line) and outer(dashed line) galactic disks (i.e at 
z — 0) for a uniform source distribution in the galactic disk are shown in Figs. 
5&6 together with direct measurements in the solar vicinity. Proton data points are 
taken from [31-38] and the electron data are from [39-45]. The injection rates of 
the particles from the sources are such that all the model spectra are normalized to 
the observed spectra at lOGeV since the solar modulation effects are believed to 
be negligible for particle energies > lOGeV. For protons, it can be seen that the 
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spectrum in the inner and outer regions almost overlap each other when normalised 
to observations at lOGeV. But, for electrons the spectrum in the inner region is 
found to show a slight steepening compared to that in the outer region. The values 
of the model parameters used in the calculations are as follows: 

The total galactic magnetic field B is found to vary with the galactic radius R [48]. 
Its value is about 10 fiG at R = 3 kpc, and decreases to about 6 fxG near R = 8.5 
kpc, and then down to 4 fxG at around R = 17 kpc. For the present analysis, B is 
taken as 10(6) /iG for inner(outer) galaxy. For the mean density of ionised gas in the 
galaxy n HII = n HII (R)dR/(R 2 - where n HII (R) = J^n HII (R, z)dz/h, 
the density distribution nnn(R, z) given in [49] is adopted. R 2 ) takes the same 
value as in section 3.1 for calculating fin- This gives uhii = 0.05 cm~ 3 in the 
inner galaxy and n HII = 0.03 cmT 3 in the outer region. n H is taken as same as in 
section 3.1. Though the density of 2.7K cosmic microwave background radiation 
is almost constant in the universe with a value of w M br = 0-25eV/cm 3 , other 
ISRFs densities like the diffuse NIR/optical radiation and the diffuse FIR radiation 
fields do not remain constant in the galaxy, wnir varies from ~ 0.5 eV/cm 3 at 
R = 8.5 kpc [50] to ~ 2.5 eV/cm 3 at around R = 1 kpc in the galaxy. w FIR 
varies in the range of ~ 0.2- 0.3 eV/cm 3 [51] in the galactic plane. Here, a value of 
wnir = 2.0(0.5)eVy cm 3 for inner(outer) region and w FIR = 0.2eV/ cm 3 for both 
the regions are adopted. Also, the SN explosion energy is taken as E S n = 10 51 
ergs, ejecta mass as M e j= 8M@, injected particle density and energy in the SN 
shock as n = 10 -3 cm -3 and 1 KeV respectively, and an upper energy cut-off of 
WOTeV in the source term. 

It is found from the calculations that the average proton and electron spectra in both 
the inner and outer galaxy over-predict the observed data both at lower and higher 
energies. The over-prediction of the model spectra above the data in the lower en- 
ergy region may be attributed to the solar modulation effects, but that in the higher 
energy region (E> 10 GeV) where the solar modulation effects are less, clearly 
shows that both the protons and electrons spectra in the galactic disk are flatter 
than the observed spectra. This deviation of the calculated spectrum from the ob- 
served spectrum may be understood as follows: for example, because of the small 
energy-loss timescale (ti oss ) for high energy electrons compared to diffusive escape 
timescale, they cannot travel a long distance in the galaxy and only the electrons 
from nearby sources may contribute to the local spectrum. For diffusion coefficient 
D ~ (1 - 10) x 10 28 cm 2 s~\ the range of a 100 GeV electron (t loss ps 10 6 yr) is es- 
timated as R fa 2^/Dti oss (200 — 600)pc. That means only those sources within 
a distance d ps (200 — Q00)pc from the Sun can contribute to the local electron 
spectrum at 100 GeV. Detailed calculations done by Atoyan et al.1995 [29] have 
shown that the total CR electron flux near the Sun even at energies 1 — 100 GeV 
where the particles do not suffer significant energy losses, calculated assuming a 
uniform and continuous source distribution in the galactic disk in the framework 
of diffusive (without convection) propagation model, is dominated by particles in- 
jected from sources with distances d < lkpc. Therefore, it is quiet necessary to 
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know the electron spectrum generated by local sources together with the local CR 
propagation parameters, in addition to the galactic averaged spectrum (here consid- 
ered as the spectrum generated by distant sources) if the observed spectrum is to be 
understood properly. Detailed discussions about the local CR spectra (both protons 
and electrons) will be discussed in section 5. 

It should also be mentioned at this point of discussion that in the case of a station- 
ary injection of CRs into the galaxy from sources continuously distributed in the 
galactic disk, the equilibrium CR (proton and electron) spectra in the disk depends 
weakly on the actual geometry of the source distribution as long as the CR energy 
densities are normalized to the same value. Another important point is that away 
from the galactic plane, the shape of the CR electron spectrum shows a significant 
variation from that in the disk particularly in the high and low energy regions due to 
their high radiative and ionization energy losses at these energies i.e. the spectrum 
is softer by an index of Ar ~ 0.6 at a distance of lkpc away from the galactic 
plane for energies greater than ~ lOOGeV. But, the proton spectral shape is almost 
the same as in the disk except in the very low energy region due to their high ion- 
ization energy losses operating at these energies. So, in the following section while 
calculating the diffuse gamma-ray spectrum averaged over a galactic region, such 
a variation of the CR spectra with distance away from the disk will be taken care. 



4 Diffuse gamma-ray spectra 

There are three important processes for the production of diffuse non-thermal high 
energy gamma-rays in the galaxy due to the interaction of galactic CRs with the 
ISM and ISRF. They are the inverse compton (IC) and bremsstrahlung of electrons 
and the decay of 7r°-mesons produced in inelastic collisions of CR protons and 
nuclei with ambient nucleons. 



4.1 Gamma-ray emissivities 

The 7r°-mesons decay 7-rays emissivity is given by [52] and references therein as 

tf(*,Kr) = 2r. ?j Z ' E jdE„ (24) 

where E 1 ^} n {E 1 ) = E 1 + rn^ c 4: /AE 1 and m^o is the 7r°-meson rest mass energy. 
The emissivity of it mesons is calculated as 

q n o(z,E n o) = Annnr] [ dE p J p (z, E p ) d<7 ^ ' Ep (25) 

J E rnin/ E ) F rf-ELo 
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where J p (z,E p ) « (c/4n)N p (z, E p ) is the CR proton flux, da^(E^o,E p )/dE n o 
is the differential crossection for a CR proton of kinetic energy E p to produce a 
7T° with total energy E^o, rj = 1.5 is the nuclear enhancement factor which takes 
into account the contribution of heavier nuclei in CRs and ISM to the the 7-ray 
emission and E p nm (E 7T o) is the minimum proton kinetic energy that contributes to 
the production of a pion with energy E n o . 

The emissivity of IC 7-rays due to the interactions of CR electrons and the ISRFs 
using the full Klein-Nishina cross-section is given by [53] 



(z,E J ) = 8ir 2 r 2 (m 



2\ °^ph 



dE. 



J e (z,E e ) 



Eph JE th 6 El 



2q In q + (1 + 2q) (1 - q) + \ ^ 



(26) 



where 



r = 



AE ph E e 
m\ 



1 = 



Efh = — 
2 



E e T(l - E y /E e ) 
£ 7 + 



\ 



J e (z,E e 



47T 



£2 + ^X 

E p h 
N e (z,E e ) 



n p h = 



Wph 



E 



•ph 



Wph and n ph are the energy and number densities of the target photon respectively, 
Eph is the energy of the target photon, E th is the minimum electron energy that can 
scatter a target photon of energy E p h to energy E 7 , r = 2.82 x 1CT 13 cm is the 
classical electron radius and m e is the electron rest mass energy. 

The emissivity of bremsstrahlung 7-rays due to the interactions of CR electrons 
and the ISM is given by [53] 



q b ; em (z,E^)=A7Tn H / J e (z,E e ) 



dEL 



dE P . 



(27) 



where 



d<7brem(E e , -E-y) 
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Fig. 7. Averaged diffuse 7-ray spectrum in the inner galaxy (330° < I < 30°, \b\ < 5°) cal- 
culated using the CR proton and electron spectra in the region, assuming w p = lAeV/cm 3 , 
w e = 0.07eV/cm s , uh = lcm~ 3 , wmbr = 0-25eV/cm 3 , wnir = 2.0eV/cm 3 , 
WFIR = 0.2eV/cm 3 , B = WfiG and r d = 15kpc. COMPTEL and EGRET data points are 
as in [9], Milagro [56], Whipple [58], HEGRA [59,60] and TIBET [61]. Model components 
: 7T° decay (thin solid line), bremsstrahlung (dashed line), IC (dotted line), extragalactic dif- 
fuse 7- ray background (dot-dashed line) taken from [55] and total (thick solid line). 




Fig. 8. Averaged diffuse 7-ray spectrum in the outer galaxy (30° < I < 330°, \b\ < 5°) cal- 
culated using the CR proton and electron spectra in the region, assuming w p = l.OeV/ cm 3 , 
w e = 0.05eV/cm 3 , nn = 0.6cm~ 3 , wnir = 0.5eV/cm 3 , B = 6fiG, = 9.5kpc and 
all the other model parameters are same as in Fig.7. Data points: EGRET [9], Milagro [56], 
HEGRA [60] and TIBET [61]. 



E e ,min = max\E^,E{\ 

where a = is the fine structure constant, z is the atomic number of the target 
atom and E L is the lowest energy in the electron spectrum. 
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4.2 Total gamma-ray flux 



Using the 7-ray emissivities listed above in section 4.1, the total 7-ray flux with 
energy E 1 in a given direction denoted by longitude (/) and latitude (b) is calculated 

as 

1,(1, b, E,) = ^J [q*°(z, E,) + qfiz, E,) + q b ^ m (z, E,)]dr (29) 

The integration is over the line of sight distance r measured from the Sun in the 
direction of (/, b) and z = r sin b denotes the height above the galactic plane. 
Notice that 1,(1, b, E,) is independent of / in the present model since both the CR 
proton and electron spectra are assumed to be constant in the galactic disk (z = 0) 
as discussed in section 3. For calculating the it decay and bremsstrahlung 7-rays 
fluxes, the mean ISM hydrogen atom density is taken as n H = lcm" 3 in the in- 
ner galaxy and n H = 0.6cm~ 3 in the outer region (as discussed in section 3.1). 
For IC 7-ray flux, the target photons considered here are the 2.7 K cosmic mi- 
crowave background radiation, diffuse NIR/optical radiation and the diffuse FIR 
radiation fields. Their energy densities are taken as wmbr = 0.25 eV/cm 3 and 
wfir = 0.2 eV/cm 3 in both the galactic regions, and wnir = 2.0(0.5) eV/cm 3 in 
the inner(outer) region (see section 3.2). The overall normalization in calculating 
the 7r°-decay, bremsstrahlung and IC 7-rays fluxes goes as fin xr^x w p , fin xraxw e 
and w p h x x w e respectively where w p , w e and w p h = wmbr + wnir + wfir 
are the protons, electrons and the total target photons energy densities respectively 
and r d is the line of sight depth of the gamma-ray emission region. The total CR 
energy density of protons, helium and heavier nuclei in which the main fraction 
are protons, varies from the locally measured value of ~ 1.0 eV cm -3 at sunspot 
minimum (~ 0.8 eV cm" 3 at sunspot maximum) to a value of ~ 1.6 eV cm" 3 at 
~ 60 AU [54]. Also, the CR electron energy density has been found to vary from 
the local value of 0.05 eV cm" 3 to 0.2 eV cm" 3 at ~ 60 AU [54]. Here, for cal- 
culating the diffuse 7- ray flux in both the inner and outer regions of the galaxy, 
the particle injection rates from the sources are chosen such that the resulting 
CR proton energy density w p = 1. 4(1. 0)eV cm~ 3 and the electron energy density 
w e = 0.07(0.05)eV r cm" 3 for the inner(outer) galaxy. For the velocity integrated 
CO intensities given in [25], the molecular H 2 column density N(H 2 ) in the direc- 
tion of the inner galaxy varies approximately in the range of (0.95—3.4) x 10 22 cm" 2 
and that in the outer galaxy as ~ (0.28 — 1.33) x 10 22 cm" 2 for a CO-to-H 2 con- 
version factor X = 1.9 x 10 20 cm~ 2 (_ft'A;ms" 1 )" 1 . Here, we assume a value of 
N(H 2 ) = 1.5 x 10 22 cm" 2 for the inner galaxy and N(H 2 ) = 0.4 x 10 22 cm" 2 for 
the outer galaxy. Therefore, for the assumed molecular H 2 averaged number den- 
sity of n H , 2 = 0.33(0. 14)cm" 3 for the inner(outer) galaxy, the mean line of sight 
depths of the emission region ~ 15(9.5)A;pc in the inner(outer) galaxy. It should 
be noted that the individual normalizations fiu-, w p ,w e , w p h and may be slightly 
varied from the above mentioned values but the overall normalization should give 
the same value. 
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The diffuse 7-ray spectrum calculated using the above set of model parameters, av- 
eraged over the latitude range |6| < 5° are shown in Figs.7&8 for the inner and outer 
galaxy respectively together with the available experimental data. The EGRET data 
shown in the plots are taken from [9,55]. The Milagro point in the outer galaxy is 
obtained from [56] and for the inner galaxy, a conversion factor of 2.43 ± 0.55 from 
(40° < / < 100°) to (-40° < I < 40°) is adopted [57]. Other data points are upper 
limits set by Whipple [58], HEGRA [59,60] and TIBET [61] collaborations. In the 
inner galactic region, the model diffuse gamma-ray spectra obtained in this paper 
show a good fit to the observed data at almost all the energies except at ~ 85GeV 
where the model predicts a much higher flux than the observed data. Though there 
is a possibility that the CR spectra derived here may not represent the true spectra 
at around 85GeV, it should also be kept in mind that the effective energy range 
of EGRET is only about 30MeV — 30GeV [1] where the present model shows 
a good fit to the data. In the outer galaxy also the model spectra nicely explains 
the observed gamma-ray data, but the results obtained in this region are somewhat 
less convincing because of the limitation of significant number of data points. The 
EGRET points extend only upto 7GeV [9] in this region and beyond that we have 
only the MILAGRO point at 3.5TeV for (40° < I < 100°, \b\ < 5°) [56] together 
with the upper limits given by ground-based experiments like TIBET and HEGRA. 

According to Hunter et al. 1997 [1], the spectrum of diffuse gamma-ray emission 
in the EGRET energy range does not vary significantly with longitude or latitude 
for |6| < 10°. The overall results obtained in this paper agree quiet well with those 
results since the high energy CR proton spectrum hardly varies in the galaxy in our 
model and the 7r°-decay gamma-rays dominate the total spectrum in the EGRET 
region. However, the present model is not able to explain the weak evidence (~ 3a) 
of softening of the gamma-ray spectrum in the outer galaxy compared to the inner 
galaxy at very low latitudes |6| < 2° mentioned in their paper. 



5 Local cosmic-ray spectra 

It has been seen from the results discussed in the previous sections that the averaged 
CR spectrum calculated for the inner and outer regions of the galaxy can explain 
the observed diffuse 7-ray spectrum effectively well in both the galactic regions. 
But, one serious problem which requires detail study is why the galactic averaged 
CR spectrum which fits the diffuse 7-ray data so well, is flatter than the observed 
CR spectrum. This section describes a possible explanation for the problem. 

For CR protons in the framework of a steady 3-D convection-diffusion model, it can 
be seen from Fig. 9 that for a standard diffusion coefficient D = 0.7 x 10 28 cm 2 s~ 1 , 
convection velocity Vq = 15kms~ 1 kpc~ 1 and source index T p = 2.4, most of the 
protons below ~ lOOGeV are coming from sources located within a distance of 
d ~ 1.5kpc from the Sun. This is due to the fact that the CR diffusive timescale 
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goes as (tdiff oc 1/D(E)) with energy. So, high energy CRs take a lesser time to 
reach us than the low energy CRs. Those protons with energy less than lOOGeV 
may be either convected away from the galactic plane or may suffer significant 
energy losses due to ionization and coulomb interactions with the ISM (if their 
energies are less than ~ lOOMeV) or due to adiabatic cooling (if their energies are 
greater than ~ lOOMeV) before reaching us from their sources. So, protons below 
~ lOOGeV are not able to reach us effectively from distant sources. The present 
analysis assumes a constant diffusion coefficient throughout the galaxy as well as a 
continuous and uniform source distribution in the galactic disk injecting CR protons 
with a power-law source index T p = 2.4. But, this may always not be the case. If 
the diffusion coeffient in the local region is different from the galactic averaged 
value, then the scenario may be much different from what is shown in Fig. 9. One 
such scenario is shown in Fig. 10 where the diffusion coefficient is taken as D = 
10 27 cm 2 s" 1 for region within 1.5kpc from us and D = 5 x 10 28 cm 2 s _1 for region 
beyond l.bkpc. The total flux (solid line) is splitted into the contributions from local 
sources within l.hkpc (dashed line) and distant sources beyond l.hkpc (dotted line). 
It shows that the majority of the CR protons upto ~ lOTeV are liberated mainly 
by sources located within 1.5kpc provided the diffusion coefficient in the local 
region is allowed to differ appreciably from the galactic averaged value. In fact, 
studies assuming a localized nature of the galactic CRs requires the value of the 
diffusion coefficient to be D ~ (10 26 — I0 27 )cm 2 s~ 1 considering the observed CR 
anisotropy [62]. So, from now on, it will be assumed that such a variation in the 
diffusion coefficient exists in the galaxy and so the observed CRs are contributed 
mainly by local sources. 

For CR electrons, as already mentioned in section 3.4, considering the diffusive 
escape and energy loss timescales, the majority of the locally observed electrons 
are mainly generated by sources located within a distance of d < 1 kpc from the 
Sun (Atoyan et al.1995 [29]). Though the authors had discussed the problem in 
the framework of a purely diffusive model without convection, the result is almost 
same in the convection-diffusion model as well. So, without going into much de- 
tail, the observed CR electrons will be asumed to be contributed mainly by sources 
located within 1.5kpc as in the case of protons. It can be noted that for any ac- 
ceptable value of D < I0 29 cm 2 s~ 1 , the sources of the observed CR electrons 
are mainly localized. It is because of the fact that in addition to the diffusive es- 
cape timescale, various other timescales like the ionization, bremsstrahlung, inverse 
compton, synchrotron, adiabatic cooling and convective timescales govern the CR 
electron spectrum. 

Now, to model the observed CR proton and electron spectra considering only the 
sources located within 1.5kpc from us, we will proceed with the usual steps fol- 
lowed in this paper as described in sections 2&3, i.e. finding the SNR age distri- 
bution (here for the local region), then use it to calculate the source spectrum and 
then propagate the CRs in the galaxy. From a compilation of several parameters 
of 231 all-type SNRs in the galaxy given in Table 1 of Jian et al. (2005) [12], it 
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Fig. 9. Decomposed CR proton spectrum near the solar region in the framework of a con- 
vection-diffusion model for Do = 0.7 x 10 28 cm 2 s~ l , a = 0.6, Vo = Ihkms^kpc' 1 
and source spectral index T p = 2.4. The solid curve (marked as 0) is the total flux decom- 
posed into the contributions from sources located at distances d > (0.5, 1.0, 1.5, 2.0)kpc 
shown by the curves marked as (0.5, 1, 1.5, 2) respectively. The total spectrum is normal- 
ized to the observed flux at WGeV and the calculation assumes nn = 1.1 lcm -3 and 
nun = 0.022cm -3 . Data points are same as in Fig.5. 
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Fig. 10. Unnormalized CR proton spectrum calculated using a convection-diffusion 
model assuming diffusion coefficient D$ = 10 27 cm 2 s _1 for the local region (i.e 
d < 1.5kpc), Dq = 5 x 10 28 cm 2 s~ 1 for the region beyond l.5kpc and convection ve- 
locity Vo = lbkms' 1 kpc^ 1 . Here, a = 0.6 and T p = 2.4. It can be seen from the plot 
that the total spectrum (solid line) upto ~ WTeV is contributed mainly by local sources 
(dashed line) located within l.bkpc. The contribution from the distant sources located be- 
yond 1.5kpc is represented by the dotted line. 

is found that there are only 10 S-type SNRs which are located within a distance 
d < l.bkpc from the Sun. Out of these, 8 SNRs have known ages. The /, b, d and 



23 



Table 2 

Parameters of 8 S-type SNRs with d < 1.5kpc : 



I 

(°) 


b 

(°) 


d 

(kpc) 


Age(t) 
(yrs) 


logiot 


65.3 


5.7 


1.0 


14000 


4.146 


73.9 


0.9 


1.3 


10000 


4.000 


74.0 


-8.5 


0.4 


14000 


4.146 


89.0 


4.7 


0.8 


19000 


4.279 


114.3 


0.3 


0.7 


41000 


4.613 


119.5 


10.2 


1.4 


24500 


4.389 


160.9 


2.6 


1.0 


7700 


3.886 


299.2 


-2.9 


0.5 


5000 


3.699 



age(t) parameters of these 8 SNRs are listed down in Table 2. Column 5 lists the 
log w t of each SNR. A gaussian fitting on the distribution of these SNRs ages using 
equation (2) gives the fitted parameters as : a = 0.415 ±0.236, ai = 4. 153 ±0.124 
and a 2 = 0.278 ± 0.142. Using this age distribution for the nearby SNRs, the CR 
source spectrum is obtained using equation (6) and is found to be T p = 2.24 ± 0.12 
for protons with E » M p and T e = 2.27 ± 0.08 for electrons. The large error bar 
on the spectral indices are due to the small sample size (only 8 SNRs with known 
ages) involve in the analysis. After propagating the CRs in the galaxy, the proton 
and electron fluxes at the solar position can be obtained using equations (16) and 
(21) respectively by taking r x = and r 2 = l.hkpc. The energy spectrum thus cal- 
culated are shown in Figs. 11&12 for CR protons and electrons respectively. The 
calculation assumes a continuous and uniform source distribution, a diffusion coef- 
ficient D = 0.3 x 10 27 cm 2 s" 1 and a convection velocity of V = l^kms^kpc' 1 
in the local galactic region. The solid line represents the model spectrum calcu- 
lated for a sharp energy cut-off at lOOTeV. For electrons, the model spectra calcu- 
lated for an exponential energy cut-off at 5TeV (dotted line) and at lOTeV^ (dashed 
line) are also shown. All the model spectra are normalized to the observed flux at 
lOGeV. The over-prediction of the model spectra above the observed spectra be- 
low E < 10 GeV can be due to the presence of solar modulation effects. In the 
case of protons, the slight deficit of the model spectrum from the observed spec- 
trum for energies greater than ~ lOTeV^ shows the required contribution of the 
distant sources as expected as discussed above (see Fig. 10) in this section, for the 
allowed variation of the local diffusion coefficient from the typical galactic value. 
The model parameters used in the calculations are chosen at (R, z) = (8.5, 0)kpc 
as fiu = nni + 2nH 2 — 1-11 cm~ 3 where nni is taken from [24] and nn 2 from 
[25], n HII = 0.022 cm" 3 [49], w M br = 0.25 eV cm" 3 , w NIR = 0.5 eV cm' 3 [50], 
wfir = 0.2 eV cm~ 3 [51], SN explosion energy E S n = 10 51 ergs and ejecta mass 
M ej = SM @ . 
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Fig. 11. CR proton spectrum calculated in the convection-diffusion model by considering 
only the sources located within l.5kpc from the Sun assuming Dq = 0.3 x 10 27 cm 2 s~ l , 
a = 0.6, Vo = Ihkms^kpc' 1 , fin = 1.11 cm -3 and uhii = 0.022 cm~ 3 . The source 
spectrum is calculated as described in section 5 and takes a power-law spectral index of 
approximately T p = 2.24 ± 0.12 for E » M p with a sharp energy cut-off at WOTeV. 
The solid line is the calculated model spectrum normalized to observations at lOGeV. Data 
points are same as in Fig.5. The slight deficit of the model spectrum from the observed 
spectrum for E > lOTeV shows the required contribution of the distant sources in the 
high energy region (see section 5). 

6 Results and discussions 



The model presented in this paper considers SNRs as the main sources of galactic 
CRs. Considering the SNR age distribution in the galaxy, the model involves the 
calculation of averaged CR source spectrum both in the inner and outer regions of 
the galaxy independent of the diffuse 7-ray data. Earlier attempts to explain the 
u GeV excess" either assumed an adhoc source spectrum optimized on the 7-ray 
data or a spectrum obtained from radio observations of SNRs which, as already 
mentioned, does not have an SNR age dependence (which actually should have 
according to shock accleration theory). After knowing the source spectra, CRs are 
allowed to propagate in the galaxy using a steady 3-D convection-diffusion model 
assuming a uniform and continuous source distribution in the galactic disk. It is 
found that the resulting CR spectra is able to fit the observed diffuse 7-ray spectrum 
quiet well both in the inner and outer regions of the galaxy. 

The present model assumes both the leptonic and nucleonic origin of diffuse gamma- 
rays in the galaxy like the models discussed in most of the earlier works. But, get- 
ting into detail one can notice that the results obtained here differ from the earlier 
results in some respects. In the inner galaxy at energies around lOOMeV — lOOGeV, 
both the 7r°-decay gamma-ray component and the IC component are important to 
explain the diffuse gamma-ray spectrum and above ~ SOOGeV^ the 7r°-decay com- 
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Fig. 12. CR electron spectrum calculated in the framework of a convection-diffusion model 
by considering only the sources located within l.5kpc from the Sun. The calculated source 
spectrum adopted here has a power-law index of approximately T e = 2.27 ± 0.08 as de- 
scribed in section 5. All the model parameters are same as in Fig. 11. The solid line is the 
calculated model spectrum for sources with a sharp energy cut-off at lOOTeV, dotted and 
dashed lines are for sources with an exponential energy cut-off at 5TeV and WTeV re- 
spectively. All the model spectra are normalized to observations at WGeV. Data points are 
the same as in Fig. 6. 

ponent mainly dominates the total spectrum, unlike in earlier models where enough 
freedom is given to the 7r° component such that it alone can possibly explain the 
GeV energy region and the IC component dominates the TeV region. In the present 
model, the normalization to the 7r°-decay component is mainly constrained by the 
flux in the TeV region whereas in the earlier models, it is mainly constrained by 
the flux in 100 MeV — lOGeV region. Here, fitting the gamma-ray spectrum in the 
range of 100 MeV — lOGeV by 7r°-decay component of radiation alone will predict 
a higher flux than the observed data in the TeV range and so the possibility of ex- 
plaining the u GeV excess" by the 7r°-decay gamma-rays alone is ruled out unless 
the proton spectrum is too steep with (T p > 2.6) which again is undesirable since 
it gives a higher gamma flux in the lOOMeV — lGeV region. In the present model, 
the underprediction of the model spectrum below the observed gamma-ray flux in 
lOOMeV — lOGeV energy region when the n° component is constrained by the 
observed flux at 3.5TeV is compensated mainly by the IC component. Also, one 
very important result that can be concluded is that the present model demands the 
CR proton spectrum to have an exponential cut-off at ~ lOOTeV so that the model 
spectrum does not exceed the observed flux in the TeV region. This value of proton 
cut-off energy agrees well with the predictions of shock acceleration in supernova 
remnants [e.g. 10]. In addition to this, the studies presented here is also able to 
explain the steepening of the observed CR protons and electrons spectra from the 
galactic averaged spectra provided the diffusion coefficient D in the local region 
is significantly less than the galactic averaged value. In short, the important results 
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of this paper can be listed as 



(i) In the inner region of the galaxy, the averaged CR proton source spectrum can 
be represented by a broken power-law in kinetic energy (E) with index T p ~ 
2.12±0.01 for E » M p and 1.66±0.01 for E « M p and that of electron by 
a power-law index T e ~ 2.126±0.004 for E » M e . 

(ii) In the outer region of the galaxy, the source spectral indices for proton goes as 
approximately 2.1 1±0.01 for E » M p and 1.66±0.01 for E « M p , and that of 
electron as 2.124±0.004 for E » M e . 

(iii) The averaged diffuse gamma-ray spectrum calculated using the galactic CR 
spectra explains the observed diffuse 7-ray spectrum in both the inner and outer 
regions of the galaxy quiet well. However, the weak evidence (~ 3a) of softening 
of the gamma-ray spectrum in the outer galaxy compared to the inner galaxy at 
\b\ < 2° [56] is still unexplained by the present model. 

(iv) The possible explaination of u GeV excess" by the 7r°-decay gamma-rays alone 
particularly in the inner galaxy is ruled out unless the proton spectrum is too steep 
with (T p > 2.6) which again is undesirable since it predicts a higher gamma flux 
in the lOOMeV - lGeV region. 

(v) Using a gaussian function fitted to the age distribution of only 8 S-type SNRs 
with known ages located within a distance of l.bkpc from the Sun, the source spec- 
tra in the local region are calculated as T p ~ 2. 24 ±0.1 2 for protons with E » M p 
and T e ~ 2.27 ± 0.08 for CR electrons. 

(vi) These local source spectra easily explain the steepening of the observed CR 
proton and electron spectra from the galactic averaged spectra provided the diffu- 
sion coefficient D is much less than the galactic averaged value. 

(vii) The present studies support the SNR origin of galactic CRs. 

Though the model presented here is able to explain the u GeV excess" problem 
quiet well and at the same time maintaining a local CR spectrum which fits the ob- 
served CR data, the results given in this paper involve some uncertainty. Firstly, the 
errors involve in the estimation of SNR ages (known ages) are quiet large and only 
for a few SNRs, the ages have a good accuracy. Secondly, in the present analysis 
the contribution to the galactic CRs from Plerion (P-) type and Composite (C-) type 
SNRs (which make a total of ~ 23% of total galactic SNRs) are not considered. 
Lastly, the large error on the local CR source spectral indices are due to small sam- 
ple size in the local region (here only 8 SNRs) and can be reduced only by using 
a statistically large sample. In fact, the number of SNRs detected is quiet less than 
the total number of SNRs expected. The total number of SNRs in the galaxy can 
be estimated as N = vT where v is the galactic supernova explosion rate and T 
is the average lifetime of an SNR. For v = 1/30 yrs -1 and an assumed range of 
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lifetime T = 10 4 — 10 5 yrs, N ~ (333 — 3333). This number strongly depends on 
the assumed lifetime and for the average value, T = 5 x 10 4 yrs, iV ~ 1667. But, 
a total of only 231 SNRs are given in Green's catalog [11], i.e. only ~ 14% of the 
expected total number of SNRs have been detected in radio. So, roughly speaking 
a total of ~ 60 SNRs are expected within d < 1.5kpc if the SNRs are assumed to 
be uniformly distributed in the galaxy. 

In future work, studies considering the contribution of P- type and C- type SNRs 
and also the diffuse emission in the low energy region will be addressed in detail. 
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